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Abstract 

The analytical approach developed by us for the calculation of the phase diagram 
for the Anderson localization via disorder [J.Phys.: Condens. Matter 14, 13777 
(2002)] is generalized here to the case of a strong magnetic field when q subbands 
{q = 1,2,3) arise. It is shown that in a line with the generally accepted point of 
view, each subband is characterized by a critical point with a divergent localization 
length ^ which reveals anomaly in energy and disorder parameters. These critical 
points belong to the phase coexistence area which cannot be interpreted by means of 
numerical investigations. The reason for this is a logical incompleteness of the algo- 
rithm used for analysis of a computer modelling for finite systems in the parameter 
range where the finite-size scaling is no longer valid. 
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1 Introduction 



In the series of our recent papers [T][2][3llHl5] an exact analytic solution for the 
Anderson localization [6] was presented. As it has been commented [71l8ll9] . 
most of our results contradict the generally-accepted rules of the scaling the- 
ory and numerical modelling. Two particular results were highly unexpected 
in the Anderson community: (i) existence of the Anderson transition in two 
dimensional (2-D) disordered systems of noninteracting electrons, and (ii) the 
Anderson transition in N-D (A^ > 2) dimensions is of first-order, and localized 
and conducting states could co-exist. 



Email address: kuzovkov@latiiet.lv (V.N. Kuzovkov). 



Preprint submitted to Elsevier 



Received 11 June 2008 



The result (i) formally does not contradict experiment [10], where the unex- 
pected presence of a metallic phase in 2-D was observed [11]. It was noted 
flOj that experiments of the last decade do not support the prevailing point of 
view that there can be no metallic state or metal-insulating transition in a 2-D 
system. The physics behind these observations at present is not understood. 

The problem is that the experimental facts present not a strong argument 
for many theoreticans who believe that they are able to calculate disordered 
systems of noninteracting electrons [T^ . This problem was studied in the liter- 
ature by various methods: scaling theory of localization, perturbation theory, 
numerical modelling, etc. The conclusion has been drawn (but not proved) 
that a phase transition in experimental systems is related most likely to the 
interactions (Coulomb, spin-orbit) neglected in theory. However, this state- 
ment is valid only for theoretical Hamiltonian-based methods. In particular, 
the scaling theory of localization [13] is a phenomenological approach and thus 
could potentially be able to take into account strong particle interactions. The 
fact that this theory [13] is unable to describe experimental data indicates its 
serious problems. It is also possible that a 2-D metal-insulator transition is 
also possible for noninteracting electrons [T|5] : however, existing theoretical 
methods are not suited for solving such a complex problem. 

It should be also noted here that the above-mentioned consensus that in 2-D 
all states are localized should be reconsidered in the light of the conflicting 
situation. Indeed, as it was noted in a recent review article [12] devoted to 
the numerical investigations of the Anderson localization, practically all nu- 
merical results are in contradiction with the analytical predictions (see also 
discussion in Ref. |14j). Despite author [12] claims that results of numerical 
analysis could be treated as obtained "from first principles", one should re- 
member that any numerical analysis is based on certain analytical predictions. 
In particular, the necessity to modify the numerical algorithm for finite-size 
scaling was discussed [S] ■ The applicability of finite-size scaling studies has also 
been questioned [T51IT6] . In other words, there is an obvious conflict between 
the results of different analytical and numerical methods. 

This is why instead of the discussing point (i) [existence vs nonexistence of 2-D 
phase transition] it is more reasonable to consider consequences of point (ii) in 
the light of the applicability of finite-size scaling. This is important since this 
is a key method in both numerical investigations of the Anderson localization 
|12] and analytical realization of the finite-size scaling [8j. The concept of a 
finite-size scaling is taken from theory of phase transitions [T7lll8|19j . and in 
fact is based on the extrapolation of the results for finite-size systems (length 
L < oo) for the thermodynamic limit, L = oo. It is clear that as the result 
of such an extrapolation, a single limiting value should be obtained. However, 
this is true only for the second-order phase transitions and in fact, the finite- 
size scaling [T7lll8fl9j could be applied only to this particular case, but not to 
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the theory of phase transitions in general. 



This leads usually to the assumption of the existence of the correlation length 
which is divergent at the critical point [12j. The method cannot be used for the 
first-order transitions with a known multiplicity of solutions p1[3]|5] where for 
an arbitrary set of system parameters two phases can coexist (localized and 
delocalized in our case). As mentioned in Ref. [1], a trivial example is the phase 
co-existence of ice and water. Only properties of the two pure homogeneous 
phases are physically meaningful, whereas an extrapolation of the results for 
a heterophase system has no real meaning. However, when the result of finite- 
size scaling do not fit to the second-order transition, the quick (and wrong) 
conclusion is often drawn that no phase transition occurs at all [1] neglecting 
a possibility of the first-order transition. 

The type of the phase transition does not follow from any first-principles [20] 
and should not be the matter of speculations [H] . This is defined by the exact 
solution (if it exists) or careful analysis of information taking into account 
all possible alternatives. From the experimental point of view, the first-order 
phase transitions are observed relatively easily. In particular, the direct electro- 
static probing [21] and photoluminescence spectroscopy [22] clearly show co- 
existence of localized and metallic regions associated with 2-D metal-insulating 
transition. New theories are put forward to address this issue |23ll2l] . 



In the present paper, we apply our approach [T][2][S1IHIS] to a 2-D system in an 
external magnetic field. The field could be included into a 2-D Hamiltonian by 
means of a Peierls factor [12]. This model is used for numerical analysis of the 
critical Hall regime in disordered systems in a strong perpendicular magnetic 
field. The prevailing viewpoint for this model is that no metal-insulator tran- 
sition can occur here. It is believed also that there is the critical energy Ec in 
each Landau band at which electrons are delocalized [25] . In other words, this 
is nothing but a generalization of our dispute on phase transitions [TP|5|^f5] 
for the case of an external magnetic field. We demonstrate below that the 
model can be exactly solved using the analytical approach [T|2][3] for a cer- 
tain range of magnitude of the control parameter (magnetic flux). As a result, 
the generally-accepted viewpoint on the metal-insulator transitions should be 
revised. Moreover, unlike the standard opinion [2S] that the delocalization at 
the critical energy Ec arises entirely due to a strong magnetic field, we show 
below that the existence of phase transitions without magnetic field and the 
critical energy E^ in the magnetic field are mutually related. 



Note that under the exact solution in this paper we mean only the calculation 
of the phase diagram with the localization length. To solve this problem, we 
use the recursive equation for the Cauchy problem with fixed initial condi- 
tions. There are some limitations of the analytical theory. An exact solution 
is only possible for the conventional Anderson model with a diagonal disor- 
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der, where on-site potentials are independently and identically distributed. 
We do not calculate here transport and other (energy spectrum, etc) proper- 
ties, which corresponds to the problem with fixed boundary conditions (the 
Dirichlet problem) and has no exact solution. 

The paper is organized as follows. In Section 2 the Schrodinger equation is 
discussed for noninteracting particles in the presence of a perpendicular mag- 
netic field. The calculation of the joint correlators is discussed which permits 
to extract analytically information on the Lyapunov exponents and to draw 
the phase diagram. Since we extend our approach P^EfS] for the case of mag- 
netic field, we discuss mostly the necessary modifications of our approach. This 
is why knowledge of our previous papers is prerequisite here. The Lyapunov 
exponents are found using theory of functions of complex variables (via search 
of poles of the H{z) function on the complex plane). Two conformal mappings 
are discussed in Section 3 which permit us to reduce the pole search to a rela- 
tively simple algebraic problem. In Section 4 the main results are summarized 
and nontrivial aspects of our theory analyzed. 



2 Main definitions 

2.1 The model 

We consider a single-band disordered Anderson tight-binding model with Schro- 
dinger equation [12] 

Ipl+l^m + 1pl-l,m + exp(-ifco0^i,m+l + exp(iA;oO^«,m-l = {E - ei^rn)i'l,m{l) 

describing the properties of noninteracting particles in the presence of a per- 
pendicular magnetic field. The magnetic field enters the transfer terms con- 
necting nearest neighbors via the phases, where = B{eay/h is the magnetic 
fiux through the elementary plaquette of the size (the lattice constant a 
and the hopping matrix element are equal to unity). The on-site potentials ei^m 
are independently and identically distributed with existing first two moments, 
(^i,m) = and {ei^rn^i'^i) = (^"^61^116,^,1, where the parameter a characterizes the 
disorder level. 

All moments of the disorder distribution define some physical properties (en- 
ergy spectrum, etc), however, phase diagram under study is defined only be 
the second moment, this is an exact result p!|2|3] . In this paper, along with 
the phase diagram, we determine also the Lyapunov exponent 7, which is the 
inverse of the localization length. The Lyapunov exponent is a typical order 
parameter. It is well known in the theory of phase transitions that there is 
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no unambiguous definition of an order parameter. Tliis also liolds in our case. 
Many different definitions are possible p!|2|3] . also dependent on all moments. 
We use here the ip"^ -definition [Tll^|[5] . where the Lyapunov exponent depends 
only on the second moment. This is why calculation of the two first moments 
of the disorder distribution is sufficient for solving our problem. 

Without disorder {ei^m = 0), the energy spectrum can be obtained analytically 
for ko = 27Tp/q, where p and q are coprime integers (rational number of flux 
quanta per plaquette). In system with an external magnetic field each energy 
band will split into several Landau subbands (g subbands). Unlike 3-D case, 
these subbands are not overlapped with each other 

In a study of disordered system {ei^m 7^ 0) we discuss a general method, 
whereas in calculations restrict our self to the cases q = 1,2,3, assuming p = 1. 
For q = 1 exp(ifco/) = 1, and the relevant Schrodinger equation coincides 
with that without the magnetic field (only one Landau subband). That is, 
the problem is reduced to the results [Ij indicating the existence of a 2-D 
metal-insulator transition. The peculiarity of the case of two Landau subbands 
(g = 2) is a lack of the subband gap in the energy spectrum of the ordered 
system. For g = 3 three subbands are separated by finite width gaps. There is 
no qualitative difference between g = 3 and larger g values. 

The peculiarities of the model ([T]) caused by the magnetic field are: (a) the 
equation contains in general case complex coefficients and one has to look for 
complex solutions, (b) the two directions in a 2-D system (characterized by 
indeces / and m) formally enter Eq. ([1]) asymmetrically which could be inter- 
preted as anisotropy, (3) the coefficients exp(±2/co0 entering the Schrodinger 
equation are periodic functions. This is why the magnetic field problem solu- 
tion requires a generalization of the mathematical formalism p^|l2|l3|B] . to be 
discussed below. 

2.2 Signals and filter function 

In our previous papers [Tf^f^ the Anderson localization was considered as the 
particular case of the generalized diffusion. This process, as any diffusion, is 
characterized by the divergence of the average over an ensemble. For instance, 
the localization effects lead to the divergence of the second momentum (di- 
agonal correlators), {\ipn,m\'^), as the function of the index n {discrete time), 
77. — *■ 00. The divergence is characterized by the Lyapunov exponent 7, which 
is the inverse of the localization length, ^ = I/7. Hereafter the average (. . .) 
means the ensemble average over random potential realizations. 

An existence of the so-called fundamental mode responsible for the correlator 
(|^n,mP) divergence was strictly proved in [3]. It was shown that the Anderson 
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localization problem could be reduced in the general case to the so-called 
signal theory characterizing the filter function hn (localization operator). The 
complementary parametrs to hn are input signals s^°^ and output signals s„ 
[L2.3j . These quantities fulfil the fundamental equation 



n 
1=0 



(0) 



(2) 



The input signal characterizes the ideal system with a zero disorder parameter 
a. Its basic properties are: it is bounded, Is^^-*! < oo, as n — >• oo, for physical 
(delocalized) band solutions, and is divergent for the formal (non-physical) 
solutions. For physical s^^ solutions, Eq. ([2]) describes the transformation of 
the input signal into another signal s„, called output. The latter charac- 
terizes the disordered system with a 7^ 0. If for a given energy E and disorder 
parameter a the output signal is bounded, < 00, as n — >• 00, this indi- 
cates transformation of input physical solutions into delocalized ones. In the 
opposite case, |sn| — > 00, as n — > 00, the output corresponds to the localized 
states. 

It was shown [T |[2|3] that the divergence of the output signal with a nonzero 
disorder is independent of the input signal, but is a fundamental feature of the 
localization operator hn (the filter function or system function). As n — 00, 
\hn\ < 00 for delocalized states, but it is divergent for localized solutions, 
\hn\ 00, n 00. 

Use of the Z-transform 

Hi^) = £ ^ (3) 

n=0 ^ 



(similar for both input and output signals) permits to reduce the problem to 
a search of the poles of the function H{z) of the complex argument z: 

H-\z) = 0. (4) 



The mo 

dified5(°)(^), S{z) andH{ z) remain to be called the input and output 
signals, and the filter function, respectively, whereas Eq. ([2]) reduces 

S{z) = H{z)S^''\z). (5) 



It was shown p!|2|3] that knowledge of the fundamental characteristics of the 
system - the filter H{z) - permits to determine the phase diagram of the sys- 
tem (regions of localized and delocalized solutions), as well as the localization 
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length. On the contrary, the input and output signals by themselves are not 
of great interest. However, the filter, input and output signals are complemen- 
tary characteristics: to find the filter H{z), one has to solve analytically exact 
equations for signals. This is possible under certain conditions, but the solu- 
tion is very complicated [3j- However, the calculations could be considerably 
simplified after a careful analysis of the structure of the solution [3]. 



2.3 The block structure of the Schrddinger equation 
Let us assume periodicity of the coefficients 

exp(z/cog) = 1, (6) 

and introduce two integer indices: n = 0, 1, ... (block number or discrete time) 
and j = 0,l,...,g — 1 (coordinate in the block), provided I = nq + j. 

Let us define new amplitudes as 

<^n,L = i^nq+j,m, (7) 

and the operators C^^^ which acts on the index m 



£(i)^(i) ^ ^^(i) _ 6"'''°^'"^' 09^^ 



ikojm' ^^(j) 

m+m' 
m'=±l 



Due to periodicity, = C^'^ . 

Since the coefficients in Eq.([T]) are complex variables, in order to calculate 
correlators, we introduce the complex conjugated functions Tp'£>^ and operators 



v=±i 

Let us rewrite Eq. ([1]) in a recursion form. Inside block (j = 1, 2, . . . , g — 2) 

one gets 

rn,m n,mrn,m ' rn,m rn,m ' V "7 

whereas on its boundaries (j = and j = q — 1) 
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,(1) =_p(o) ^(0) ,Ao) (0) nr 



,m ' 



Similar relations are found for complex-conjugated functions. In particular, an 
analog of Eq. (fTOj) reads 



wr^=-^i^s+^'V3-^&^^- (13) 

2.4 Correlators 

Following our approach |i3j, the diagonal correlators (j^n]m\'^') should be calcu- 
lated. To this end, a complete set of equations for the two types of correlators 
has to be analyzed: 

^^^■n^W = (^gL^S>' (14) 
y^'\nu, = U^lTfi!;'^) . (15) 



Here (^J,/^ = <^^^_il, x^^\n)^^i = x^^\n)i^rn- 

The equations for correlators can be derived relatively easily. Say, to define 
x^''Kn)m,i one has to multiply the left and right sides of the above-given equa- 
tions, in particular, Eqs. ffTOj) and f|T3|) are used inside the same block. Then, 
the obtained relations should be averaged over the ensemble of random po- 
tentials: 



x^^^'\nU, = (egL^i) (v'gL^S) + (16) 



n,m 



When doing so, the causality principle piPfS] is taken into account, which 
means that all amplitudes v^^l^, Tp^n ] on the r.h.s. of obtained equation are sta- 
tistically independent with respect to the potentials e^^l^, eg]. In other words, 
when calculating the average quantities on the r.h.s. of equations, only ampli- 
tude correlations are essential. The calculation of averages over the potentials 
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is quite trivial: we use know averages for the first two moments, (^Sn]^^ = 
and 

{^n^m^n'!n,') = 0^^'^j,i"5n,n'^m,m', (17) 



provided a remains the disorder parameter. 

Proceeding this way, one gets on the r.h.s. of the derived equations along with 
the correlators x^^\n)m,i sought for, also additional terms with the correlators 
y^^\'n)m.i- The equations for y^^\n)m,i could be derived in a similar way, r.h.s. 
are expressed through the correlators x^^\n)m,i. and y^^\n)m,i- That is, the 
equation set is complete. 

Note that the boundary conditions for the amplitude in the general form, 
ipQ,m = 'ipi^m = ttm are reduced to the boundary conditions for the correla- 
tors, e.g. a;W(0)„,z = 0, x(i)(0)„,; = amdh Due to an arbitrary choice of the 
field the translational invariance does not occur. However, the complete 
analytical solution of the correlator equations is possible by means of the dou- 
ble Fourier transform and Z-transform [3j. The same is true for the problem 
under study. This approach is exact but rather lengthy. 

The derivation could be, however, considerably reduced, taking into account 
the important result of our paper [3]. In fact, we introduced the correlators 
as the tool in our study, but all we need is only the term in the diagonal 
correlators x^^\n)m,m called the fundamental mode; it is divergent for the 
delocalized states. Namely this quantity corresponds to signals s„ in Eq. ([2]). 

The fundamental mode is invariant with respect to the argument shift in the 
boundary condition a'^ = am+mo- The equations for signals could be easily ob- 
tained by the formal replacement of the boundary conditions in equations for 
the correlators [3]: x^^\0)m,i = CimPii is replaced for x^^'>{Q)m,i = amCii = ^m-i- 
Here the function corresponds to the Fourier transform T{k) = 

After such a replacement the problem becomes translationally invariant, the 
correlators depend only on the argument divergence r = m — I 



x^^\n)m,i ^ x^^\n)r, (18) 
y^'\n)m,i ^ y^'\n)r, (19) 

whereas the diagonal correlators become independent on the m-argument (fun- 
damental mode) and transform into a set of signals 

x^'\n)m,m^ s{j,n). (20) 
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For simplicity we use hereafter the correlator symmetrization which does not 
change the diagonal elements sought for 

x^^\n)r = ^[x^'\n)r + x^^\n)^r], (21) 
y^'\n)r = l[y'^^\n)r + y'-^\n)^r]. (22) 
By definition, 

x^'\n)o = s{j,n). (23) 
2.5 The Z-tmnsform and the Fourier transform 



The analytical solution is based on the use of the two type of algebraic trans- 
forms: the Z-transform 



S(j-) = f:^. (24) 

n=0 ^ 



(^)^ ('25) 



and the Fourier transform, 



fO)(^)^= J ^X(^\z,k)e-"''. (26) 

J ZTT 



The Eq. (1231) for the correlator-signal relation transforms into 

x^^\z)o = Sij,z), (27) 

which gives also the following useful integral relation: 

J^X^^\z,k) = S{j,z). (28) 

— TV 

2. 6 Equations 

As a result, we obtain the equation set for the correlators (j = 0, 1, . . . , g — 1) 
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zX^^+^\z, k) - T{k)Sj^i = a^S{j, z) + (29) 
^z-^X^^~^\z, k) + C]{k)X'^^\z, k) - 2Cj{k)Y'^^\z, k), 

(^'+1) {z, k) = {z, k) + C,{k)X^^^ {z,k). (30) 

Here the functions 

£j{k) = E -2cos{k - jko) (31) 



arise from the Fourier transform of the operators C^^\ The periodic conditions 
take place with respect to the index j: X^'^^z, k) = X^'^\z, k). 

The solution algorithm is quite simple. For a given q a set of 2q equations 
(1291) . (!30!) has to be solved. As a result, one gets q relations for the correlators 
X^^\z, k) which are linear in signals S{j, z) and function T{k). After use of Eq. 
(1251) we obtain the set of linear equations for q signals S'(j, z). The coefficients 
here are expressed through integrals. Taking into account coefficient symme- 
try with respect to index transformation, all partial signals S{j, z) could be 
combined into the total signal 



S{z) = Y.S{3,z). (32) 
i=o 



In all cases the equation for S{z) is reduced to a general form, Eq.([5]), and 
H{z) = 1 as cr = takes place. As a result, we derive equations for the 
filter function H{z) and the input signal S^^\z). Since for calculating the 
phase diagram the filter function H{z) is sufficient, we do not discuss here 
the bulky equations for the input signals S^^\z). The more so, the analysis of 
input signals S^^\z) [3] leads to the trivial conclusions: the input signals are 
bounded if the wave functions without disorder form the band, and unbound 
outside the band. However, the existence range of the band without disorder 
could be found by means of standard methods, without complicated analysis 
of the asymptotic behaviour of the signals. 



2.7 Main relations 



Let us consider three cases, q = 1,2,3. For q = I one gets fco = 27r, Coi^k) = 
C{k), where 

£.{k) = E-cosik). (33) 
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To illustrate the derivation method of the equations for S{z), we consider a 
simple case of one Landau subband {q = 1). For other cases, only the final 
relations are presented. 

For q = 1, let us denote X'^^^z, k) = X{z, k), Y'^^^z, k) = Y{z, k), respectively. 
One gets the equation set 

{z - z~^)X{z, k) - T{k) = a^S{z) + C\k)X{z, k) - 2C{k)Y{z, k), (34) 

{z+l)Y{z,k) = C{k)X{z,k). (35) 

This set is solved trivially with respect to X{z, k)\ 



(^ + 1) 



C\k) 



X{z, k) = T{k) + a'^Siz). 



The equation for the signal S{z) could be derived using the Eq. 



(36) 



^X{z,k) = S{z). 



(37) 



The resulting equations were discussed in our Ref.[T]. 

The results for q = 1, 2, 3 could be presented in a general form 



1 



1-a' 



zi + l) } dk[Ni{E,z) - N2{E,z)C{k)] 



H{z) ~ ' (z'i-l) J 27r {zi + 1)7^9 - £2(A;) 

— TT 

The coefficient Ni{E,z), N2{E,z) here are defined below. 
For g = 2, 3 C{k) reads: 



(3J 



C{k) = Coik)Ciik)-2, 



(39) 



Cik) = C,{k)C^{k)C^{k) - [C^{k) + C,{k) + C^{k)]. 



(40) 



Simple trigonometrical transformations show that 



£,{k) = 6-2cos{qk). 



(41) 
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The parameter S = S{E) equals S = E {q = 1), S = E"^ — 4 {q = 2) and 
S = E{E'^ — 6) {q = 3). Irrespective of the q value, the parameter 6 G [—4,4] 
for each subband in an ideal system (cr = 0). 

Notice that the integrand in Eq. (155]) contain always Ni{E, z) — N2{E, z)C{k) 
linear in C{k). This simplification arises due to our regular use of symmetry 
and trigonometrical transformations. 

In particular, for q = 2 the intermediate calculations contain the following 
intergrals 



J = / * m (42) 

— IT 

It could be easily shown that the integral does not depend on the index j. In 
particular, in the integral numerator, Eq. fH21) . Cj{k) could be replaced by 



^ ^ E = 2E'-2- C{k). (43) 

^ j=0 



Similarly, in the case q = 3, one can use the result invariance with respect to 
the cyclic index transpositions: 



Cj{k)£{k) C,{k)C{k) = Em, (44) 

^?W^^E^?(^) = 2 + ^^ (45) 
{Ci{k)C2{k) - 1)2 ^ 6 + 3E^ - UE^ - 2EC{k). (46) 



As a result, one obtains the coefficients in the numerator of the integral, Eq. 
(I3HD. For g = 1 one gets Ni{E, z) = 1, N2{E, z) = 0. For g = 2 



iVi(E, z) = i£!±ll + 2(^2 _ (47) 

z 



Lastly, in the case q = 3 
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N,{E, z) = (6 - 14^2 + 3^^) + (z^ + ^-2) + (2 + E^){z + z-^), (49) 

iV,(E,.) = 2E^:i|^. (50) 



3 Conformal mapping 

3. 1 The first conformal mapping 

Note that the input and output signals are real values. According to signal 
theory [27j, the roots Zi of Eq.(j4l) either lie on the real axis or arise as complex 
conjugate pairs, Zi,Zi. The properties of the filter function H{z) and physical 
interpretation of the mathematical solution are determined by the location of 
the largest (per modulus) root of Eq. (jl]), Zmax = max{|zj|}. As Zmax > 1, the 
filter H{z) is unstable and corresponds to the localized states, otherwise for 
Zmax < 1 {stable filter) and delocalized states occur. This is why, to determine 
the stability region boundaries of the filter H{z), we restrict ourselves to the 
z values in the upper half-plane, < arg z < n. 

Search of the complex roots of Eq. (jl]) could be simplified by performing two 
conformal mappings. 

First, we change the complex variable z for the new parameter ( = z'^. 

(1) With a magnetic field for ko = 2tt and without magnetic field ( formally, 
g = 1) this is a trivial transformation, C, = z, so that the C, variable is 
defined in the upper half-plane. 

(2) For q = 2 we use the relation ( = z"^, so the ( variable is defined on a 
whole complex plane, < arg ( < 2tt. The inverse transformation gives 
z = respectively. 

(3) In the case g = 3 we divide the upper half-plane where z is defined, 
into two regions, (a) In the first region, where < arg z < 2n/3 we 
define ( = z^, so that ( again is defined on a whole complex plane < 
arg ( < 271. The inverse transformation yields z = C^^^. (b) In the second 
region, where 27r/3 < arg z < vr we also assume ( = z^ , but now define 
the inverse transformation as z = C}^^ ■ exp(z27r/3). In this region the 
variable C, is defined in the upper half-plane. 

As a result, we obtain ^-presentation of the filter function H{(). Next we seek 
the roots of the equation H~^{() = and obtain for g = 1, 2, 3 
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1 



1 - a 



2 (C + 1) fdk [N,{E,z)-N2{E,z)Cik)] 



(51) 



3.2 Second conformal mapping 

When calculating intergrals in Eq. (ISTl) . it is convenient to use the relation 
well-known for functions of complex variables 



J 271 T ± 2 cos{qk) i^^ - ' ^ ' 

This relation holds for arbitrary g = 1, 2, . . . provided Imr > (the complex 
variable r is defined in the upper half-plane) . 

The integral in Eq. (l5T|) defines the function which is non-analytic on the unit 
circle |C| = 1 [T|2] . When performing the inverse Z-transform, this leads to a 
double solution. It is convenient [Tp] to perform a second conformal mapping 



where the choice of a sign defines one of the two filter functions, H±{w). For 



w = +{C'^^'^ + C "^^^) (or for w = —{('^^'^ + C '^^'^))^ the outer (inner) part of 



the unit circle, \(\ = 1, transforms onto the upper half-plane w = u + iv, 
V = Imw > 0. The circle itselfs maps onto the interval u G [—2, 2] on the real 
axis (f = 0). The transformation inverse to Eq. flS^ reads pQ 



3.3 Parametric representation of the pole diagram 

As shown [T]|2] . physical information on the localization could be obtained from 
the H±{w) using standard methods for the functions of complex variables. The 
key issue here is the location of poles on the complex plane (pole diagram) 




w = ±(C^/^ + C'^') 



(53) 



w"^ , wi r. ^ 

C = -1 + — ± — V4- w2. 
^ 2 2 



(54) 



m 



Let us present H±[w) in the following form 



H. 



{w) = l-a^R±{w). 



(55) 
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The main idea is quite simple. The function H±{w) has its poles where 



a^R±{w) = 1 (56) 

takes place. So, we have the parametric w-representation of the pole diagram. 
The relevant expressions for R±{w) are quite simple and allow analytical so- 
lution. In particular, for g = 2 one gets 



R+(w) 



R-{w) 



1 

1 

V4 - 



4 + 6 
{6-wy-4 

A + 6 



+ 



2 + 6 + w 
2 + 5~w 



A -{6 + wy ^ J A -{5-w) 



(57) 
(58) 



The diagrammatic technique for seeking roots of Eq. (156!) was described in 
detail in appendix of Ref. [2]. This is why we present below only the main 
results. 



4 Results 



Stable filter 



The function H^{w) (so-called stable filter [27]) defines the existence region 
of the delocalized states. General results of its study (Appendix in Ref. [2]) 
can be summarized as follows. The H^{w) allow physical interpretation if the 
equation a'^R_{w) = 1 either has no roots (which is possible only for space 
dimensions higher than two), or the roots of the equation, w = u + iv, are real 
and lie in the interval w = u (z [—2, 2]. 

That is, in our case it is sufficient to find the existence region of the real roots 
w of the equation a^R_{w) = 1 in the interval w = u E [—2, 2], provided the 
absence of other (complex) roots. As a result, one can determine the function 
a = aQ{6) (or a = cro(-E), when changing the 6 variable for E) which gives 
the critical disorder destroying the localized states. This function is plotted 
below, Fig.l, when analyzing the pase diagrams for g = 1, 2, 3. 

For the illustration we present some results for the case q = 2 (two Landau 
subbands) which were obtained partly analytically, partly using numerical 
methods. The equation a'^R-{w) = 1 for 6 = has two real roots in the 
interval w = u E [—2,2] provided the disorder parameter < a < ao[5 = 
0) = 0.828. As cr > ao{6 = 0), equation has complex roots having no physical 
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interpretation. As already was mentioned, the interval S G [—4, 4] corresponds 
to the old band (for zero disorder). In the intervals —4 < S < 6q {6q = —0.741) 
and 2 < 5 < 4 the function o"o((5) = 0, that is infinitesimally weak disorder 
already destroys all delocalized states. However, the delocalized states are 
allowed in the interval < 5 < 2, provided the disorder parameter a is below 
some critical value. Note, that the function (To(5) is multi- value in the region 
5 e ((5o, 0) and has a characteristic triangle shape. The delocalized states are 
destroyed in the beginning even due to very weak disorder, as the a parameter 
grows, formation of the delocalized state phase is possible. The point 5q on 
the energy scale {5 = E'^ — 4) corresponds to the \E\ = Eq = 1.805. The 
formation of a phase of delocalized states is possible in the energy interval 
\E\ G (1.805,2.449), respectively. 



4-2 Unstable filter 



The localization is characterized by the unstable filter Hj^{w) [2]. The main 
result is quite simple [T|2] : we are seeking one specific root of the equation 
a'^R^{w) = 1. This root corresponds to the pole maximally distant from the 
coordinate origin, zq = Zmax, and could be related to the localization length 
^ via a simple equation: zo = exp(27) [1] (where the Lyapunov exponent 7 = 
1/^). After two conformal mappings the point zq maps onto wq = 2cosh(2g7) 
and the mathematical problem is reduced to seeking a real root of the equation 
a'^R+{w) = 1. 

It could be shown that the root wq sought lies on the real axis in the interval 
w = u> uq where 



uo = 2 + \5\. (59) 

The value of cr = corresponds to the wq = uq, and wq monotonically in- 
creases, as a increases. That is, in the limiting case of a — > the Lyapunov 
exponent 7 = 7(5, a) becomes 

7(5,0) = lcosh-i(l + |5|/2). (60) 
2q 



The minimal value of Mq = 2 occurs for S = 0, that is the critical energy values 
Ec {one per each subband) being defined by the equation 6{E) = are well 
defined {anomaly discussed below). Here ■j{S = 0,0)= 0, i.e. the localization 
length ^ is divergent at the critical points Ec, which are equal to: Ec = 
{q = 1), E, = ±y/2 {q = 2), and E, = 0, ±^6 {q = 3). 
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In the limit of a small a one gets 



7(5 = 0, a) oc a. (61) 

From Eq. fl60|) for \5\ <^ 1 one gets 

7(5,0)oc.\^, (62) 

Omitting the proof details, let us discuss the results for the Lyapunov exponent 
7(-E, a) (i.e. using E instead of parameter 6) as shown in Fig.2, for g = 1, 2, 3. 

Eq. (!62|) suggests that in the vicinity of the critical point Ef. (taking into 
account \5\ oc \E — Ec\) in the limit a ^ one gets 7(-E, 0) oc \E — E^^^'^. That 
is, we expect a strong anomaly in the Lyapunov exponent '~f{E, a) behaviour. 
In surface plot, Fig.2, this anomaly looks like deep canyons around the points 
[E = Ec,cr = 0). The canyon banks are very steep, therefore minimal values 
of the Lyapunov exponents lie on the canyon bottom, i.e. form the line. This 
anomaly completely confirms the general point of view that the localization 
length diverges at a single energy point at the center of the Landau band [25j . 



4-3 Phase diagram 



As was mentioned, a study of the if_ (w)-function allows to determine the 
existence region for the delocalized states, whereas H^{w) characterizes the 
localized states. It was found [2] , however, that the localized states exist at all 
energies and non-zero disorder parameters. 



localized states 




localized states 



phase co-existence^ 




localized states 



phase co-existence^ 




Fig. 1. Phase diagram for g = 1 (a), g = 2 (b) and q = 3 (c). 



Therefore, in order to draw the phase diagram the analysis of the stable filter 
H^{w), Fig.l, is sufficient. The {unlimited) existence region of the localized 
states overlaps with that for the delocalized states {limited in both energy and 
disorder parameter). A number of such region corresponds to a number q of 
the Landau subbands, see Fig.l. Each region is located near a center of old 
Landau subbands which existed without disorder. 



18 




Fig. 2. Surface plot ^{E,a) ioi q = 1 (a), q = 2 (h) and q = 3 (c). 



As a result, the phase diagram contains one-phase regions (only localized 
states) and two-phase co-existence regions. The delocalized state phase is de- 
fined by zero Lyapunov exponent which needs no graphical illustrations. On 
the contrary, the localized state phase is characterized by the function 7(-E, a) 
shown in Fig. 2 in form of the surface plot. Note also that the results in Fig.l 
and Fig. 2 correspond to the thermodynamic limit. 

The generally-accepted point of view [25] is that in the presence of a strong 
perpendicular magnetic field the extended states appear in the centers of 
disordered-broadened Landau bands and give rise to the integer quantum Hall 
effect. Instead, we obtained a more complex situation. In each Landau sub- 
band there is a region of phase co-existence. This region contains the critical 
points {E = Ec, cr = 0) where the localization length is divergent. Moreover, 
the observed anomalies (the canyon bottom in the surface plot of 7(-E, a) in 
Fig. 2) also lie in the region of phase co-existence. 

Since the traditional viewpoint is based mostly on the results of numerical 
investigations, it is unclear how justified is here the thermodynamic limit 
transformation. It was reasonably noted [25] that we have to keep in mind that 
numerical calculations always deal with finite systems, while, strictly speaking, 
phase transitions exist only in infinite systems. It was also believed there that 
finite-size scaling theory could serve as theoretical framework to analyze data 
for finite systems, in order to extrapolate results to infinite system size. Our 
analysis of the phase diagram allows to understand better this confiicting 
situation. 

First of all, the general result should be mentioned. Both Fig.l and Fig.2 show 
that despite a quantitative difference, as a number of the Landau subbands 
q increases, the same qualitative picture remains as it is observed for q = 1. 
There is the co-existence region in each subband including the critical point 
Ec and the relevant anomaly in the Lyapunov exponent (localization length). 
However, as q increases, the region size systematically decreases. One can ex- 
pect that in systems with large q the two-phase region practically disappears. 
This is why there is some reasoning behind the statement [T2] that the model 
([1]) does not exhibit a metal-insulator transition. This is true indeed but only 
for a large number of Landau subbands. 
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The case of g = 1 considered above allows us to establish also a link be- 
tween systems with and without external magnetic field. As was noted, under 
a strong magnetic field but for p = q = 1, Eq.([T]) coincides formally with that 
without magnetic field. Therefore Fig. la and Fig. 2a illustrate also the phase 
diagram and the 7(-E', cr)-surface plot without magnetic field. In other words, 
our correction of the statement p^|l25] that in the presence of a strong mag- 
netic field extended states appear in the centers of Landau subbands, namely, 
a demonstration that the extended states could arise in fact in a wider region 
of the phase-coexistence, sheds also light on theoretical interpretation of the 
delocalized state existence without the magnetic field. Our important conclu- 
sion is that an emergence of the critical points is not a new effect caused by 
the magnetic field but the fundamental feature of a system even without the 
imposed field. 

In the main region denoted in Fig.l as the "localized states" the situation 
is simple and clear: a single phase exists for a given E and a parameters. 
Since the analytical results correspond to the average over the random po- 
tential realizations, the existence of the delocalized states is here impossible. 
The random potential realizations could give slightly different results but this 
means, however, nothing but homophase fluctuations. All quantities are well 
deflned and the terminology of the self- averaging quantities is justifled here. 
This region of the phase diagram is not problematic also for standard numeri- 
cal investigations. Indeed, the numerical study deals with the fluctuations in a 
homophase flnite-size system, with the main focus on extrapolation of these re- 
sults with typically bad statistics to the thermodynamic limit. Despite the fact 
that flnite-size scaling is not strictly justifled, this faces no fundamental objec- 
tions. The problem is that the flnite-size scaling works well in the cases when 
the region of the phase diagram under study contains points with a divergent 
correlation length. These points correspond to the relation ^{Ec, a = 0) = 0. 
However, as follows from Fig. 2, such points are absent in the homophase re- 
gion. The more so, the values of 7(-E, o") are not small since they lie outside 
the canyon bottom. 

When considering the phase co-existence region in the phase diagram, theory 
of the flrst-order phase transition is commonly used. It is obvious that use of 
the two parameters, E and a, does not deflne uniquely the phase state of a sys- 
tem: an ensemble of random potential ei^m fluctuations could be heterophase, 
not homophase ones. A study of such ensembles faces a serious theoretical 
problem. It is common in experimental studies to characterize a large sample 
in terms of a single phase, giving no statistical analysis. However, in theoreti- 
cal studies the macroscopic features are associated with the ensemble average. 
It is obvious that the averaging over the heterophase ensemble has no physical 
sense: this is true only for characteristics of individual homogeneous phases, 
self- averaging quantities lose here any sense pLj. 
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The standard idea that an increase of the sample length L brings us to the 
situation where the self-averaging effect automatically takes care of the en- 
semble statistical fluctuations and the statistical error can be estimated from 
sampling different realizations of the disorder [2S] does not work here. An ad- 
ditional problem is that our approach p!|2][3] suggest no method to estimate 
the weight (probability) of individual phases. 

Therefore, the phase co-existence region contains the critical points {E = 
Ec, a = 0) where the localization length is divergent. However, the finite-size 
scaling cannot be used here, since the situation differs qualitatively from that 
for the second-order phase transition. That is, observation of the critical points 
does not necessarily justify use of the finite-size scaling. 

The very fact of the heterophase ensemble existence produces fundamental 
difficulties for numerical investigations. In finite systems the phase concept is 
not defined, thus this is impossible to separated statistical contributions of 
the two phases. Indeed, how could one distinguish large but rare homophase 
fluctuation on a given parameter in a finite (typical small) size and the het- 
erogeneous fluctuation, especially when the statistics is quite poor (relatively 
small number of random potential realizations)? We are not able to suggest 
any new algorithm for analysis of the relevant numerical investigations, we 
can only discuss the consequences of the heterogeneous ensemble existence. 



4-4 Homophase interpretation 

The possibility of phase co-existence is neglected tradicionally in numerical 
investigations [12]. It is believed that even if the phase diagram consists of 
the two regions with two phases, in each of these regions the relation between 
parameters E, a and the phase is unique. In other words, an existence of only 
homophase ensembles is assumed, which could arise in the case of the second- 
order phase transitions as well as without transitions. Let us discuss what 
could be a result of the substitution of a heterophase ensemble for homophase 
one, while performing numerical investigations. We will demonstrate that the 
general statement follows uniquely from the homophase interpretation. 

1) In the phase co-existence region both realizations with delocalized state, 
7^0, and localized state, 7 > 0, exist. As was said above these contributions 
cannot be separated (definition of a phase exists only for infinite systems [25j), 
and the statistical analysis leads to a formally calculated average (7) > 0. 
This is equivalent to the statement that the system does not exhibit a metal- 
insulator transition (formally (7) > means existence of the localized states 
only). Therefore the homophase interpretation of 2-D numerical investigations 
with and without magnetic field leads to well-known conclusions [T^PS] 
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2) At the same time, the heterophase fluctuations should manifest themselves 
in the statistics as abnormally strong fluctuations of the 7 parameter, which is 
unusual for homophase systems and bring into question the treatment of the 
7 parameter as the self-averaging quantity. The situation is complicated by 
the presence in the same region of the critical points {E = Ec, cr = 0) where 
the localization length diverges. Since such critical points arise typically in the 
second-order transitions, the quick conclusion suggests this interpretation with 
strong homophase fluctuations nearly the critical points instead of the real 
heterophase fluctuations. Moreveor, the observation of critical points could be 
also used to justify the finite-size scaling. 

3) The anomaly of the Lyapunov exponent 7(-E, cr) near the critical points is 
well pronounced in Fig. 2. This specific anomaly determined by the numerical 
methods, with a well pronounced peculiarity of the canyon bottom in the plot 
'y{E,a) permits a traditional interpretation [25] about critical points and the 
delocahzed states. 

The analysis above clearly demonstrates that the thermodynamic limit has 
no correct solution by means of numerical investigation. The extrapolation in 
terms of the finite-size scaling assumes an existence of the homophase ensem- 
ble. However, there is no proof of the the existence of such an ensemble within 
the same method. The heterophase idea is not constructive since it suggests 
no alternative to the finite-size scaling. The importance of exact analytical 
solutions in such conflicting situations is self-evident. 



5 Conclusion 

We have shown that the contradiction between the results of the present ana- 
lytical approach (see also previous papers p!|2|3Pf5] ) with those of numerical 
investigations [T21I25] arises due to different interpretations of the phase di- 
agram. Our result [T1[^IIB| correspond to the thermodynamic limit where the 
phase concept is well defined [20] • It was also shown that the Anderson lo- 
calization problem is characterized by the phase diagram with specific phase 
co-existence regions which need analysis in terms of first-order transitions. In 
those regions where the ensemble of heterophase fluctuations occurs the idea 
of self-averaging quantities fails. Analysis in the phase co-existence region is 
impossible for the finite systems since the phase concept is not defined here. 
Therefore, our analytical results [T][2][3llHl5] cannot be reproduced by means of 
numerical investigations since we cannot suggest any new algorithms for the 
numerical result analysis. 

In its turn, the main numerical investigation results p^1[25] were obtained 
namely for the finite-size systems. Transition to the thermodynamic limit is 
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traditionally performed in numerical investigations by means of the finite-size 
scaling. The latter assumes (directly and indirectly) existence of the ensemble 
of homophase fluctuations which permits result extrapolation to the infinite 
single-phase system. This approach fails in the case of two phase co-existence. 
Namely a wrong use of the hypothesis of a single phase in the heterophase 
case leads to incorrect theoretical conclusion about phase-transition absence in 
2-D systems which obviously contradicts experimental data p^lITT] . Simultane- 
ously, even in the framework of incorrect homophase interpretation numerical 
methods are able to detect phase co-existence regions, interpreting these as 
the critical points with relevant delocalized states. 
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